Vanishing quasiparticle density in a hybrid Al/Cu/Al single-electron transistor 
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The achievable fidelity of many nanoelectronic devices based on superconducting aluminum is lim- 
ited by either the density of residual nonequilibrium quasiparticles n qp or the density of quasiparticle 
states in the gap, characterized by Dynes parameter 7. We infer upper bounds n qp < 0.033 fim~ 3 
and 7 < 1.6 x 10 -7 from transport measurements performed on Al/Cu/Al single-electron transistors, 
improving previous results by an order of magnitude. Owing to efficient microwave shielding and 
quasiparticle relaxation, typical number of quasiparticles in the superconducting leads is zero. 



I. INTRODUCTION 

Active research and debate on the origin and density of 
residual quasiparticles in aluminum-based superconduct- 
ing quantum circuits prevails currently^. Aluminum 
is a widely used metal in the field of low-temperature 
mesoscale electronics due to its superconducting prop- 
erties and its tendency to form a native oxide that can 
be employed as a tunnel barrier. The practical perfor- 
mance of superconducting devices is often degraded by 
excess quasiparticle processes that do not follow from 
the assumption of full thermal equilibrium and Bardeen- 
Cooper-Schrieffer (BCS) form for the quasiparticle den- 
sity of states of the superconducting electrodes. The fig- 
ures of merit that we will address in this work are the den- 
sity of non-equilibrium quasiparticles n qp and the Dynes 
parameter 7 for the normalized density of quasiparticle 
states in the gap. These parameters are presently ma- 
jor limiting factors for the coherence time of Josephson 
junction qubits 2 ^ , relaxation time of highly sensitive ra- 
diation detectors^, the ultimate temperature reachable 
by normal-insulator-superconductor (NIS) junction re- 
frigerators 6 , and potentially for the accuracy of the SINIS 
turnstile-, a contender for the realization of a metrologi- 
cally accurate source of quantized electric current. 

The lowest values for 7 that have been obtained from 
subgap conductance measurements of NIS junctions and 
SINIS single-electron transistors (SETs) are 2 x 1CT 5 
and 1 x 10~ 6 , respectively^. Subgap conductance of 
defect-free NIS junctions is dominated by two-particle 
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Andreev tunneling-. For n qp , studies on supercon- 
ducting qubits 2 -^ and resonators^ have reported a low- 
temperature saturation to 10—55 /mi -3 due to an uniden- 
tified excitation mechanism. Recently, a significantly im- 
proved result n qp < 1 /im~ 3 has been obtained by a 
transmon qubit realization 9 . A similar number can be 
inferred from the quasiparticle tunneling rates reported 
in Ref. [if] for a superconducting single-electron transis- 
tor (SET) design with normal metal quasiparticle traps. 
In this work, we present quasiparticle transport measure- 
ments of Al/AlOx/Cu SETs combining carefully imple- 
mented shielding against external black body radiation 
and quasiparticle traps. The results yield unprecedented 
upper bounds n qp < 0.033 /im~ 3 and 7 < 1.6 x 10~ 7 . 
The experiment was performed on samples similar to 
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FIG. 1. (a) Large scale electron micrograph of a sample 
showing the electrical connections to the active region in the 
middle. Portions of the ground plane, covered by an insu- 
lating AlOx layer of 25 nm thickness, are also visible, (b) 
Active area of one of the measured samples, consisting of two 
Al/AlOx/Cu/AlOx/Al SETs that are coupled capacitively by 
a Cr wire located underneath the AlOx layer. Colored arrows 
illustrate the monitored in- and out-tunneling events (green 
and red, respectively), and the macroscopic electrometer cur- 
rent (yellow) in a configuration where the upper SET is used 
as the electrometer. Shadow-evaporated leads terminate at 
ohmic Au/Al contacts beginning 10 /J.m away from the junc- 
tions (not shown). 



that pictured in Fig. [TJ We monitored the charge state on 
the central island of a hybrid, i. e., SINIS-type, SET using 
another hybrid SET as an electrometer. The electrome- 
ter current idet was read out using a room-temperature 
current amplifier in the so-called DC-SET configuration. 
Most of the capacitive coupling between the SETs was 
provided through a 7 /im long Cr wire galvanically iso- 
lated from the SET metal layers. Following the technique 
introduced in Ref. 0, all electrical leads of the lxl mm 2 
pattern on the chip from the bonding pads to the active 
region were capacitively shunted by a conducting ground- 
plane that was electrically isolated from the leads by the 
25 nm AlOx layer. Pathways for quasiparticle trapping 
from the aluminum electrodes are provided by the over- 
lap through the tunnel barrier oxide to the copper elec- 
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trodes extending to within few hundreds of nm of the 
tunnel junctions, and the ohmic contact between Al and 
Au films beginning 10 /im away from the junctions. We 
also fabricated and measured a reference sample with- 
out the ohmic Al/Au contact, but having an otherwise 
equivalent design. 



the life-time broadened Dynes density of states 1 ^ for the 
quasiparticle excitation spectrum of the S electrode 



n s (E) 



Re- 



E/A + ij 



y/(E/A + i 7 ) 2 - 1 



(4) 



II. THEORETICAL DESCRIPTION OF 
SINGLE-ELECTRON PROCESSES 

The basis of our theoretical modeling of the charge 
transfer is the Golden rule expression for the first-order 
tunneling rate 
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is the quasiparticle den- 



sity of states in the superconducting electrode, fs (/n) 
is the occupation factor in the superconducting (nor- 
mal) electrode, and P(E) is the probability to emit 
energy E to the electromagnetic environment during 
the tunneling. The occupation factors are taken to be 
Fermi functions at temperature Tn (Ts), i.e., /n,s(-^) = 
[l + exp(E/k B T NiS )}~ 1 . The Golden rule formula with 
P{E) = 6(E), equivalent to a zero- impedance environ- 
ment, has been succesfully used to describe a wide range 
of charge and energy transport phenomena in SINIS 
structures^ 1 - in the range E > A. Next, we will con- 
sider theoretically the mechanisms that can cause excess 
quasiparticle processes in the subgap range E < A. 

(i) Even if charge transport is completely described 
by the above model, quasiparticle thermalization in both 
the N and S electrodes can be nontrivial at sub-kelvin 
temperatures due to strongly suppressed electron-phonon 
coupling. In a later section, we demonstrate that quasi- 
particle temperature Tn does not deviate significantly 
from the bath temperature To at temperatures above 
50 mK. On the other hand, we characterize excess quasi- 
particle excitations in the superconductor by the density 
of nonequilibrium quasiparticles 



'qp 



2D(E F ) / dEn s (E)f s (E), (2) 

J A 



where D(Ep) denotes the density of states at the 
Fermi energy. We use the literature valued D(Ep) — 
1.45 x 10 47 m _3 J _1 . At base temperature, we have 
A/^bTn) ~ 50, and hence the induced quasiparticle 
tunneling in the gap assumes a bias-independent rate 
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(ii) Experimentally observed finite subgap conduc- 
tance in NIS junctions is ofted modeled by introducing 



where parameter 7 effectively expresses the quasiparticle 
density of states in the middle of the gap as a fraction 
of the density in the normal state. Equation (j4]) may 
not be valid far from the gap edges 1 ^, and we stress that 
in the present work the Dynes model is only used as a 
tool to assess the effect of subgap states on quasiparticle 
tunneling. 

(Hi) The P(E) function appearing in Eq. (JlJ can be 
calculated from the autocorrelation function of phase 
fluctuations over the junction 1 ^. In previous worki^, it 
has been explicitly demonstrated that photon assisted 
tunneling (PAT) due to microwave irradiation (/ > 
A/ft, = 50 GHz) originating from outside the sample 
stage can be a dominant factor in the dynamics of metal- 
lic single-electron devices. In the present experiment, a 
known source of harmful phase fluctuations is the high- 
frequency component of detector back-action. It origi- 
nates from the switching noise that results from load- 
ing and unloading the detector SET island as the probe 
current is transported through. In addition, blackbody 
radiation from higher-temperature stages of the cryostat 
can reach the junction due to insufficient filtering of the 
signal lines or leaks in the radiation shields enclosing the 
sample. The approximation P(E) = ^^§r^, valid for 
E < and sufficiently weak Sy^, gives a straightfor- 
ward relation between the absorptive part of the P(E) 
and the power spectrum Sv(w) of the voltage noise over 
the junction. The above-described detector back-action 
is an instance of random telegraph noise (RTN) , for which 
the relevant high-frequency part of the power spectrum 

can be written as Sy ct (uj) — ( ^asr ) ^eS*> where k 
is the fraction of island charge coupled from the device 
under test (DUT) to the detector, C^ ot is the total ca- 
pacitance of the detector SET island, and £(w) describes 
high-frequency attenuation of the Cr wire^ 7 -. 

(iv) Certain higher-order processes could appear in 
the experimental detector traces as single-electron pro- 
cesses. The processes in question are Andreev tunnel- 
ing of a Cooper pair into the island followed by a rapid 
relaxation of a single quasiparticle through one-electron 
tunneling, and Cooper pair-electron cotunneling, which 
changes the charge on the island by one electron. The- 
oretical predictions for the rates can be made based on 
the results of Ref . [H using the known sample parameters 
and a value g/N — 10~ 5 for the normalized conductance 
per channel^. Predicted higher-order rates are orders of 
magnitude smaller than the experimental observations 
for bias voltages |Vd s | < 100 /xV. 



3 



1000 



3100 



10 - 




-200 -150 -100 -50 50 100 150 200 



FIG. 2. Bandwidth-corrected single-electron tunneling rates 
as a function of Vds at charge degeneracy measured in dif- 
ferent setups at base temperature: PT (18 mK, open up- 
ward triangles), PDR1 (50 mK, open circles), and PDR2 
(50 mK, open squares). For PDR1, data from higher tem- 
peratures is also given: 131 mK (filled stars), and 158 mK 
(filled diamonds). All these measurements were performed 
with the same sample employing Au/Al contacts for quasipar- 
ticle trapping. The two thin solid lines represent thermally 
activated rates calculated for the known sample parameters at 
158 mK and 135 mK. The dashed line is the theoretical rate 
for 7 = 1.6 x 10~ 7 at Ts = Tk = 50 mK, and the horizontal 
thick line represents the rate induced by n qp — 0.033 /tm~ 3 . 
Open downward triangles: Base temperature data from a ref- 
erence sample without Al/Au contacts measured in PDR1. 
For ease of comparison, tunneling rates from the reference 
sample have been scaled by the ratio of junction conductances 
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FIG. 3. (a) Measured width of the Coulomb staircase steps 
as a function of bath temperature for different bias voltages 
of the DUT. (b) Observed zero-bias tunneling rates at base 
temperature for different detector currents, and a linear fit 
to the data. The plotted current is the mean value of the 
detector trace from which the transition rate was determined. 



III. EXPERIMENTAL METHODS 

To determine experimentally the single-electron tun- 
neling rate for a fixed bias voltage Vd s , the gate volt- 
ages were adjusted so that the detector was at a charge 
sensitive operating point, while the DUT was at charge 
degeneracy. At degeneracy, owing to sufficient E c of 
the DUT, only the two degenerate charge states have 
a non-vanishing population. Hence, the expected elec- 
trometer signal is RTN, and the transition rate is given 
by T = r le (eVds/2) + r le (-eV ds /2), where two identical 
junctions have been assumed. The state transitions were 
identified from the recorded electrometer traces by dig- 
ital low-pass filtering followed by a threshold detector. 
We used the analytical model of Ref. to compensate 
for missed transitions due to finite detector bandwidth. 

To test the coupling of microwave radiation to the junc- 
tions, we repeated the experiments in several cryostats 
equipped with different wiring and shielding solutions. 
Two independent setups, henceforth denoted by PDR1 
and PT, to yield the lowest tunneling rates, suggesting 
that the external microwave radiation was suppressed to 
a negligible level. In both setups, standard solutions were 
used for the microwave filtering of the signal lines: 1 m 
of Thermocoax in PDR1 and a combination of Thermo- 
coax and powder filters in PT. However, similar to other 
recent works^i, the sample chip was protected against ra- 
diation shining from the higher-temperature parts of the 
cryostats by two nested rf-tight shields. The microwave 
shields were attached to the sample stage body by either 
threads sealed with indium (PDR1), or by screws (PT). 
For illustration, we present data also from setup PDR2 
that is similar to PDR1, but in which the Thermocoax 
lines terminate to a connector that was not sufficiently 
rf-tight. 



IV. RESULTS AND CONCLUSIONS 

Our main experimental data, the measured transition 
rates as a function of the device bias Vd s , are presented 
in Fig. [3J A saturation temperature below which the ob- 
served rates did not decrease was found around 80 mK. 
Using the value R T = 1.1 Mft (2.0 Mfi and 25.0 Mfi 
for the junctions of the more asymmetric reference sam- 
ple) obtained from device I-V characteristic measured 
with an ordinary ammeter, we determined the value 
A = 210 fiV by fitting the 158 mK data to the ther- 
mally activated (TA) rates assuming Tn = Tg = Tq. The 
large scale I-V characteristic is consistent with this value. 
Data obtained at 131 mK agrees with TA predictions for 
135 mK except in the range |Vd s | < 25 fiV where the 
TA rates become comparable to the saturation floor. In 
contrast, the theoretical TA rates for 50 mK are below 
0.01 Hz in the Vds range that was accessible in the count- 
ing experiment. 

In order to aid the analysis of the base temperature 
results, we present two auxiliary results ruling out plau- 
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sible explanations for the observed low-temperature sat- 
uration. First, to demonstrate the thermalization of the 
normal metal island to 50 mK, we study the width of the 
transition between two charge states as the gate charge 
of the DUT is swept past a degeneracy point, following 
the procedure of Ref. O The obtained step widths in 
PT cryostat are presented in Fig. GH^a). For bias voltages 
far below the gap edge, anomalous broadening of the step 
at low temperatures is a clear indication of non-thermal 
processes. However, at the highest studied bias voltage 
Vd s = 350 ^iV, the linear TA regime extends down to 
50 mK. Biasing of the DUT does not affect the tempera- 
ture of the normal metal island as the electronic cooling 
power for |Vd s | < 350 /iV is negligible^. 

Secondly, we studied experimentally the influence of 
detector back-action by measuring zero-bias rates at the 
base temperature at different operating points of the de- 
tector. A linear dependence of T on /dot can be observed 
in the data from PDR1 cryostat presented in Fig.^b) as 
predicted by the P(E) theory for the small detector cur- 
rents employed here (/det "C eA/h). By extrapolating to 
zero detector current, we deduce that for the typical de- 
tector current of 30 pA in the time traces on which Fig. [5] 
is based, the contribution from detector back-action is 2- 
3 Hz, which is a significant but not a dominating fraction. 

The tunneling rate data from the reference sample 
without Al/Au contacts shown in Fig. [5] displays a clear 
plateau consistent with a quasiparticle density of n qp 
ii.fi') fim~ 3 . For this sample, the dominant quasiparti- 
cle relaxation channel at base temperature is tunneling 
through the oxide barrier to the normal metal shadow. 
This allows us to infer a homogenous injection rate of 
r qp = 3 x 10 5 s /zm , presumably due to external 
radiation. Applying the quasiparticle diffusion model of 
Ref. [2 with a thermal energy distribution, the Al/Au con- 
tacts are expected to bring about a f 00-fold reduction in 
n qp at the junctions. From the observed rates in PDRf 
and PT setups, we can infer the previously stated upper 
bound n qp < 0.033 /Ltm -3 for the sample with Al/Au con- 
tacts as illustrated in the figure, i. e., a 20-fold reduction. 
Given the nominal thickness of 30 nm for the aluminum 
film, we find that the expected number of quasiparti- 
cles in the total volume of one of the aluminum leads is 
less than 0.1, i.e., the superconductors are nearly free of 
quasiparticles in a time-averaged sense. 

However, based on the absence of a flat plateau in ex- 



perimental bias dependence, non-equilibrium quasipar- 
ticles cannot account for a majority of the observed 
base temperature tunneling events. The curve for Dynes 
model with j = 1.6 x 10~ 7 is a better match to the ex- 
perimental results in the range |Vd s | < 100 /iV, more so 
for the data from PT cryostat, but the present amount 
of data does not allow one to make a definite statement 
about the presence of subgap states at 7 ~ 10~ 7 level. 

Finally, it should be noted that the base tempera- 
ture tunneling rate at zero bias is four times higher in 
the imperfectly shielded PDR2 setup compared to those 
achieved in PDR1 and PT setups. Elevated rates in 
PDR2 are caused by PAT due to stray blackbody radi- 
ation. The results from PDR1 and PT are very similar, 
and part of the discrepancy between them can be at- 
tributed to uncertainty in the detector bandwidth com- 
pensation (in the notation of Ref. 20, Tdot = 2000 Hz for 
PDR1 and 125 Hz for PT) and different probing current. 
Also, for the high base temperature tunneling rates in 
bias range 100 — 200 fj,V, we find PAT to be a plausible 
explanation. The observed rates can be reproduced by 
assuming a microwave noise spectrum with an exponen- 
tial high-frequency cutoff corresponding to an effective 
temperature of 350 mK, and a spectral density of the 
order of 1 pV/VWi at / = 50 GHz. 

In conclusion, we have demonstrated that microwave 
shielding and enhanced quasiparticle relaxation play a 
key role in achieving the highest possible performance of 
superconducting aluminum based devices. The possibil- 
ity to employ these techniques in future realizations of 
quantum information processing devices presents excit- 
ing prospects. 
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